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The  incidence  of  traumatic  brain  injury  (TBI)  among  military  personnel  is  at  its  highest  point 
in  U.S.  history.  Experimental  animal  models  of  blast  have  provided  a  wealth  of  insight  into 
blast  injury.  The  mechanisms  of  neurotrauma  caused  by  blast,  however,  are  still  under 
debate.  Specifically,  it  is  unclear  whether  the  blast  shockwave  in  the  absence  of  head 
motion  is  sufficient  to  induce  brain  trauma.  In  this  study,  the  consequences  of  blast  injury 
were  investigated  in  a  rat  model  of  primary  blastTBI.  Animals  were  exposed  to  blast  shock- 
waves  with  peak  reflected  overpressures  of  either  100  or  450  kPa  (39  and  110  kPa  incident 
pressure,  respectively)  and  subsequently  underwent  a  battery  of  behavioral  tests.  Diffusion 
tensor  imaging  (DTI),  a  promising  method  to  detect  blast  injury  in  humans,  was  performed 
on  fixed  brains  to  detect  and  visualize  the  spatial  dependence  of  blast  injury.  BlastTBI 
caused  significant  deficits  in  memory  function  as  evidenced  by  the  Morris  Water  Maze, 
but  limited  emotional  deficits  as  evidenced  by  the  Open  FieldTest  and  Elevated  Plus  Maze. 
Fractional  anisotropy,  a  metric  derived  from  DTI,  revealed  significant  brain  abnormalities 
in  blast-exposed  animals.  A  significant  relationship  between  memory  deficits  and  brain 
microstructure  was  evident  in  the  hippocampus,  consistent  with  its  role  in  memory  func¬ 
tion.  The  results  provide  fundamental  insight  into  the  neurological  consequences  of  blast 
TBI,  including  the  evolution  of  injury  during  the  sub-acute  phase  and  the  spatially  dependent 
pattern  of  injury.  The  relationship  between  memory  dysfunction  and  microstructural  brain 
abnormalities  may  provide  insight  into  the  persistent  cognitive  difficulties  experienced  by 
soldiers  exposed  to  blast  neurotrauma  and  may  be  important  to  guide  therapeutic  and 
rehabilitative  efforts. 


Keywords:  traumatic  brain  injury,  blast  neurotrauma,  memory  dysfunction,  diffusion  tensor  imaging,  magnetic 
resonance  imaging 


INTRODUCTION 

The  incidence  of  traumatic  brain  injury  (TBI)  among  military 
personnel  in  modern  combat  is  the  highest  of  any  conflict  in  U.S. 
history  (1,  2).  Up  to  20%  of  combat  veterans  met  the  criteria 
for  TBI  on  post-deployment  screening  (3),  and  exposure  to  blasts 
from  improvised  explosive  devices  contributed  to  the  unprece¬ 
dented  rate  of  mild  TBI  (mTBI).  The  overwhelming  majority  of 
events  were  categorized  as  mTBI  (4),  typically  characterized  by 
a  consistent  clinical  scenario  of  postconcussive  symptoms,  cog¬ 
nitive  dysfunction,  and  other  functional  impairments  that  follow 
a  gradual  course  of  recovery  during  the  initial  days  and  weeks 
after  injury  (5, 6).  However,  increasing  evidence  suggests  that  even 
mild  head  trauma  may  result  in  postconcussive  symptoms  and 
comorbidities  many  months  after  mTBI  (7-10).  Moreover,  even  a 
single  blast  episode  may  have  long-term  pathogenic  potential  and 
cause  enduring  neurodegeneration  (11).  Understanding  the  neu¬ 
rological  consequences  of  blast  neurotrauma  is  essential  to  enable 
prevention  and  guide  therapeutic  and  rehabilitative  efforts. 

Despite  the  considerable  number  of  experimental  investiga¬ 
tions  of  blast  TBI  using  rodent  models  reported  in  recent  years, 


the  mechanisms  underlying  blast  injury  and  the  factors  affecting 
outcomes  remain  unclear.  Whole-body  or  head-only  blast  par¬ 
adigms  have  highlighted  the  potential  for  different  mechanisms 
of  brain  injury  following  blast  exposure  (12,  13).  Likewise,  the 
effects  of  the  blast  shockwave  alone  (primary  injury)  may  be  dif¬ 
ferent  than  those  involving  rapid  head  acceleration  subsequent  to 
blast  (tertiary  injury)  and  were  recently  compared  directly  (14, 
15).  Understanding  the  time-course  of  recovery  from  blast  neu¬ 
rotrauma  is  also  vital  to  relate  injury  mechanisms  with  patient 
outcomes.  Behavioral  assessments  conducted  at  acute  (12, 16-19), 
sub-acute  (20,  21),  and  chronic  time  points  (13,  20,  21)  indicate 
prolonged  neurological  abnormalities  following  blast.  However, 
given  the  differences  in  injury  protocols  and  other  experimental 
variations  between  studies  conducted  at  a  single  timepoint  after 
injury,  assessments  performed  at  several  time  points  would  enable 
a  better  characterization  of  the  time-course  of  deficits. 

Another  important  factor  to  understand  injury  tolerance  and 
clinical  outcome  is  the  effect  of  blast  magnitude  (i.e.,  shock- 
wave  overpressure).  A  number  of  experimental  protocols  incorpo¬ 
rated  differing  overpressure  severities  to  examine  the  physiological 
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consequences  of  higher  magnitude  exposures  (22-27),  but  fewer 
studies  have  assessed  whether  higher  magnitude  overpressures  lead 
to  greater  behavioral  deficits.  This  is  important  for  the  quantifi¬ 
cation  of  injury  tolerance  and  may  also  relate  to  the  duration 
of  the  recovery  curve.  While  some  studies  demonstrated  greater 
mortality  rates  at  higher  overpressures  (18),  a  lack  of  a  consistent 
relationship  between  overpressure  and  cognitive  ability  (19,  28, 
29),  or  neuromotor  performance  (18)  has  frequently  been  noted. 
A  recent  study  by  Cernak  et  al.  demonstrated  dose  dependency  in 
Rotarod  and  Open  Field  Test  (OFT)  outcomes  for  mice  exposed 
to  whole-body  103  and  190  kPa  shockwaves  (13).  However,  a  need 
exists  to  elucidate  the  dose-response  effects  and  temporal  recov¬ 
ery  of  behavioral  outcomes  following  primary  blast  shockwave 
exposure,  given  somewhat  conflicting  outcomes  from  previous 
investigations. 

Expanded  efforts  to  screen  military  personnel  for  potential 
exposure  to  mTBI  have  addressed  the  growing  burden  of  head 
trauma  in  combat.  However,  there  is  still  no  definitive  marker 
of  mTBI,  which  has  hampered  clinical  decision-making.  Non- 
invasive  imaging  techniques,  particularly  magnetic  resonance 
imaging  (MRI),  have  considerably  advanced  the  detection  and 
understanding  of  subtle  brain  changes  following  mTBI.  Diffusion 
tensor  imaging  (DTI),  an  advanced  MRI  technique,  is  a  promis¬ 
ing  method  for  non-invasively  examining  the  effects  of  TBI.  DTI 
probes  the  molecular  motion  of  water  within  living  tissues  to 
infer  the  presence  of  microscopic  structural  abnormalities.  DTI 
has  uncovered  changes  in  veterans  exposed  to  blast  TBI  (30-32) 
and  mTBI  through  other  physical  (non-blast)  mechanisms  (33, 
34).  However,  despite  its  promise,  many  questions  remain  regard¬ 
ing  the  use  of  DTI  as  a  diagnostic  marker.  While  experimental 
models  of  blast  TBI  were  essential  in  identifying  the  physiological, 
metabolic,  and  pathological  consequences  of  blast  TBI,  DTI  has 
been  infrequently  applied  to  experimental  animal  models  of  blast 
under  well-controlled  laboratory  conditions  that  would  provide 
fundamental  insight  into  the  neurotrauma  caused  by  blast  forces. 

The  objective  of  this  work  was  to  examine  the  behavioral 
changes  from  acute  to  sub-acute  time  points,  examine  the  depen¬ 
dence  of  behavioral  outcomes  on  different  levels  of  shockwave 
overpressures,  and  quantify  structural  damage  to  brain  tissues 
using  DTI  following  exposure  to  blast.  Importantly,  the  study 
employed  a  head-only  primary  blast  injury  paradigm  to  pro¬ 
vide  fundamental  insight  into  the  mechanisms  of  blast  shockwave 
neurotrauma  and  its  associated  behavioral  consequences. 

MATERIALS  AND  METHODS 

SHOCK  TUBE  APPARATUS 

A  custom  shock  tube  (Figure  1A)  with  a  3.6-cm  inner  diameter, 
3.0-m  driven  section,  and  0.3-m  driver  section  was  used  to  cre¬ 
ate  shockwaves  with  different  overpressure  magnitudes.  A  mylar 
membrane  separated  driver  and  driven  sections.  The  driver  section 
was  pressurized  with  helium  until  membrane  rupture.  The  shock 
tube  was  previously  characterized  and  was  shown  to  produce 
accurate  and  repeatable  shockwaves  across  a  range  of  shockwave 
overpressure  magnitudes  (35).  Pressure  transducers  (Piezotronics, 
Inc.,  Depew,  NY,  USA)  were  oriented  to  record  face-on  (reflected) 
pressures  at  a  sampling  rate  of  5  MHz  at  selected  locations  and  were 
placed  immediately  adjacent  to  the  head  for  all  blast  exposures. 


Additional  characterization  of  the  blast  pressures  was  confirmed 
by  placing  sensors  at  identical  locations  in  the  absence  of  a  rat, 
and  sensors  were  oriented  to  record  either  face-on  (reflected)  or 
side-on  (incident)  pressures. 

ANIMAL  PROCEDURES 

All  animal  procedures  were  approved  by  the  Institutional  Animal 
Care  and  Use  Committee  (IACUC)  at  our  institution.  All  animals 
were  allowed  access  to  food  and  water  ad  libitum  before  and  after 
shockwave  exposure.  A  2  x  2  factorial  design  was  employed  with 
blast  magnitude  and  post-injury  duration  as  independent  vari¬ 
ables.  Animals  were  exposed  to  a  single  blast  shockwave  with  either 
a  100  or  450  kPa  reflected  peak  overpressure  (Figure  IB),  hereafter 
referred  to  the  low-blast  and  high-blast  conditions,  respectively. 
Assessments  were  conducted  at  acute  (1-4  days  post  injury)  or 
chronic  (28-31  days  post  injury)  time  points,  hereafter  referred  to 
as  the  4-  and  30-days  post-injury  (DPI)  groups,  respectively.  Each 
experimental  group  consisted  of  9-14  rats.  A  separate  cohort  of 
eight  animals  served  as  sham  controls  and  underwent  identical 
procedures  without  exposure  to  shockwave. 

On  the  day  of  injury,  Sprague-Dawley  rats  were  anesthetized 
with  4%  isoflurane  in  oxygen  and  placed  in  a  nose  cone  for  con¬ 
tinuous  delivery  of  1.5%  isoflurane.  Rats  were  placed  39  and 
17  cm  from  the  shock  tube  opening  for  the  low-  and  high-blast 
shockwave  exposure  conditions,  respectively  (Figure  1A).  Since 
the  exhaust  gasses,  or  blast  wind,  can  lead  to  considerable  head 
acceleration,  animals  were  placed  40°  and  20°  lateral  to  the  shock 
tube  axis  to  limit  blast  wind  exposure  for  the  low  and  high-blast 
conditions,  respectively.  Prior  characterization  of  the  shock  tube 
demonstrated  minimal  blast  wind  effects  at  the  locations  chosen 
for  this  study  (35).  Moreover,  the  head  was  constrained  later¬ 
ally  and  inferiorly  to  prevent  head  rotational  acceleration- induced 
injury,  and  all  shockwave  exposures  were  conducted  with  the  rat 
head  perpendicular  to  the  radial  axis  from  the  shock  tube  opening 
(Figure  2A).  A  metal  cylinder  was  also  placed  around  the  body 
to  limit  shockwave  overpressure  exposure  of  the  torso.  Prior  work 
indicated  that  the  cylinder  was  effective  at  reducing  peak  over¬ 
pressures  to  <20%  of  overpressure  magnitudes  recorded  at  the 
head.  Following  shockwave  exposure,  animals  were  removed  from 
anesthesia  and  allowed  to  breathe  freely.  Animals  were  continu¬ 
ously  monitored  until  return  of  the  righting  reflex,  returned  to 
their  cages  for  recovery  and  observation  for  at  least  15  min,  and 
periodically  monitored  for  6  h  post  injury. 

Following  completion  of  the  behavioral  testing  protocol,  ani¬ 
mals  were  euthanized  with  an  overdose  of  pentobarbital  and  per¬ 
fused  through  the  left  cardiac  ventricle  with  250  ml  of  phosphate 
buffered  saline  (PBS)  followed  by  250  ml  of  4%  paraformaldehyde 
in  PBS.  Brains  were  excised  from  the  skull  and  placed  in  fixative 
for  48  h  followed  by  PBS  for  long-term  storage. 

HIGH-SPEED  VIDE0GRAPHY 

High-speed  videography  was  used  to  record  head  motion  during 
shockwave  exposure.  The  digital  video  camera  (Integrated  Design 
Tools  Inc.,  Tallahassee,  FL,  USA)  was  placed  directly  overhead,  per¬ 
pendicular  to  the  horizontal  plane  of  the  animal,  to  monitor  head 
motion  at  6000  frames  per  second.  Physical  markings  were  placed 
at  the  center  of  the  head  (i.e.,  lambda  landmark)  and  between 
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FIGURE  1  |  Blast  tube  characterization.  A  shocktube  (A)  pressurized 
with  helium  delivered  blast  shockwaves.  Rats  were  placed  at  specified 
distances  from  the  shocktube  opening  based  on  initial  characterization 
and  were  placed  off-axis  to  avoid  blast  wind.  Pressure  sensors  depict 


the  typical  reflected  pressures  (B)  experienced  during  the  blast 
shockwave.  Blast  peak  overpressures  realized  during  actual  exposures 
to  the  rat  head  (C)  were  highly  accurate  and  reproducible  (n  =  8  for 
each). 


the  eyes  (slightly  anterior  to  the  bregma).  Video  tracking  software 
(Tracker  4.75,  Open  Source  Physics)  was  used  to  track  displace¬ 
ment  of  markings  and  derive  head  linear  and  rotational  motion 
time  histories.  Angular  position  data  were  digitally  filtered  using 
a  lowpass  Butterworth  filter  with  a  500-Hz  cutoff  frequency  fol¬ 
lowed  by  numerical  differentiation  to  derive  angular  velocity  and 
acceleration. 

MORRIS  WATER  MAZE 

The  Morris  water  maze  (MWM)  Visuo-Spatial  Learning  Para¬ 
digm  (36, 37)  quantified  post-traumatic  anterograde  amnesia  and 
spatial  learning.  The  MWM  consisted  of  a  10  cm  diameter  plat¬ 
form  submerged  1  cm  below  the  surface  of  a  25  cm  deep  pool 
of  water.  Pool  diameter  was  183  cm.  Three  testing  sessions  were 
conducted  over  three  consecutive  days.  Each  session  consisted  of 
eight  separate  trials  wherein  rats  were  initially  placed  at  each  of  the 
four  cardinal  locations  (N,  E,  S,  W).  During  each  trial,  rats  were 
allowed  to  swim  until  finding  and  mounting  the  hidden  platform 
or  until  60  s  had  passed.  The  location  platform  was  randomized 
between  sessions  but  was  consistent  within  the  same  session,  and 
the  first  trial  of  each  session  was  excluded  from  analysis.  Visual 
cues,  including  those  internal  to  the  maze  and  external  cues  in  the 
room,  were  maintained  between  sessions.  A  computerized  tracking 
system  and  software  (Ethovision  V8.0,  Noldus  Information  Tech¬ 
nology,  Wageningen,  The  Netherlands)  recorded  movement  of  the 
rats  within  the  maze.  Latency  to  find  the  hidden  platform  and 


number  of  unsuccessful  trials  were  computed  from  these  record¬ 
ings.  Spatial  learning  deficits  manifest  as  greater  latencies  and  a 
higher  number  of  unsuccessful  trials. 

ELEVATED  PLUS  MAZE 

The  elevated  plus  maze  (EPM)  assessment  is  an  ethologically  rel¬ 
evant  test  that  quantified  anxiety-related  behaviors.  The  maze 
consisted  of  four  perpendicular  10  cm  x  90  cm  arms  connected 
by  a  10  cm  x  10  cm  central  platform.  One  pair  of  opposing  arms 
was  enclosed  by  40-cm  high  walls,  while  the  other  two  arms  and  the 
center  platform  were  uncovered.  The  entire  apparatus  was  com¬ 
posed  of  black  plexiglass  and  was  raised  50  cm  above  the  floor. 
Rats  were  initially  placed  on  the  central  platform  facing  one  of  the 
open  arms  and  allowed  to  explore  the  maze  for  5  min.  A  camera 
mounted  above  the  maze  and  Ethovision  software  automatically 
quantified  the  number  of  arm  changes,  time  spent  in  the  open 
arms,  and  the  number  of  head  dips.  Behaviors  associated  with 
increased  anxiety  included  decreased  time  in  the  open  arms  and 
fewer  head  dips. 

OPEN  FIELD  TEST 

The  OFT  was  used  to  examine  anxiety  and  spontaneous 
exploratory  behavior.  The  apparatus  consisted  of  a  flat,  100-cm 
diameter  circular  arena  with  20-cm  walls.  Illumination  was  pro¬ 
vided  by  overhead  lights.  Rats  were  initially  placed  in  the  center 
and  allowed  to  explore  the  arena  for  5  min.  The  center  of  the  arena 
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FIGURE  2  |  Head  kinematics  during  blast  Blast  was  directed  to  the  left 
side  of  all  animals  (A)  Landmarks  of  ink  placed  between  the  eyes  (red),  and 
on  the  head  (blue)  were  tracked  with  high-speed  video  to  obtain  head 
position  (drawn  to  scale)  during  exposure  to  the  blast  wave,  with  an 
example  of  the  high-blast  shown.  Angular  position  (B)  was  derived  by 
intersection  of  a  line  passing  through  the  two  points  compared  to  the  initial 
position  of  the  head,  and  demonstrated  only  limited  head  motion. 


consisted  of  a  66-cm  diameter  circle,  with  the  remainder  consid¬ 
ered  the  border  region.  The  distinction  was  indicated  within  the 
software  and  did  not  include  any  physical  markings.  An  overhead 
mounted  camera  along  with  Ethovision  software  automatically 
quantified  measures  including  total  distance  traveled  and  total 
time  spent  in  the  center  of  the  maze.  Behaviors  associated  with 
increased  anxiety  include  a  decreased  time  spent  in  the  center. 

STATISTICAL  ANALYSIS 

A  two-factor  Analysis  of  Variance  (ANOVA)  was  used  to  identify 
statistically  significant  interactions  for  behavioral  metrics  between 
blast  severity  and  post-injury  time  for  groups  exposed  to  blast, 
followed  by  post  hoc  tests  when  warranted.  In  the  absence  of  a 
significant  interaction,  a  one-way  ANOVA  was  conducted  with 
the  sham  group  included  to  test  for  a  significant  main  effect  of 
either  blast  magnitude  or  time  after  injury  using  separate  mod¬ 
els.  Post  hoc  comparisons  were  performed  if  the  overall  test  was 
significant.  For  all  tests,  significance  was  set  at  a  p-value  of  <0.05. 

MAGNETIC  RESONANCE  IMAGING 

A  9.4  T  horizontal  bore  Bruker  BioSpec  was  employed  for 
MRI  procedures.  Fixed  brains  were  immersed  in  susceptibility¬ 
matching  fluid  (Fomblin;  Solvay  Solexis,  NJ,  USA)  and  placed 
in  a  custom  20-mm  diameter  inductively  coupled  loop-gap  coil. 


Room  temperature  was  maintained  throughout  the  experiment. 
A  series  of  gradient  echo  scout  images  were  acquired  to  ensure 
reproducible  slice  orientations  and  placement  within  the  magnet. 
A  spin-echo  sequence  (TR  =  2000  ms,  TE  =  21  ms)  with  Stejskal- 
Tanner  pulsed  gradients  was  used  to  perform  diffusion-weighting 
along  12  non-orthogonal  directions  at  a  fo-value  of  1200  s/mm2 
and  two  non-diffusion  weighted  images  (b  =  0  s/mm2 ) .  A  double¬ 
echo  spin-echo  readout  (echo  spacing  =  3.3  ms)  was  used  to 
improve  signal  to  noise  ratio  (SNR)  by  magnitude  averaging  of 
the  individual  echo  images.  Twenty  axial  slices  at  a  thickness  of 
0.5  mm  and  an  in-plane  resolution  of  0.2  mm2  covered  the  entire 
brain  and  brainstem.  The  entire  MRI  experiment  was  completed 
in  approximately  1  h. 

MRI  DATA  ANALYSIS 

Diffusion  weighted  images  were  corrected  for  eddy  current  distor¬ 
tions,  and  the  diffusion  tensor  was  calculated  on  voxel-by-voxel 
basis  using  fMRI  Software  Library  (FSL).  For  spatial  registration,  a 
study-specific  template  was  created  using  tensor-based  registration 
implemented  in  the  DTI-TK  software  package  (38).  An  iterative 
procedure  using  rigid  body,  affine,  and  diffeomorphic  registra¬ 
tion  was  used  to  create  a  template  at  a  final  resolution  of  200  |xm3. 
Tensor  volumes  were  spatially  smoothed  with  an  anisotropic  three- 
dimensional  filter  implemented  in  Matlab  (39).  Maps  of  fractional 
anisotropy  (FA)  were  derived  and  used  for  subsequent  statistical 
analysis. 

Differences  between  blast-exposed  groups  and  the  sham  group 
were  assessed  using  a  two-sample  f-test.  Permutation-based  infer¬ 
ence  testing  implemented  in  the  FSL  was  used  for  statistical 
hypothesis  testing,  and  voxels  with  p  <  0.05  (uncorrected  for 
multiple  comparisons)  and  clusters  >100  voxels  were  considered 
statistically  significant.  Relationships  between  FA  and  behavioral 
tests  were  examined  using  linear  regression  analysis  of  the  blast- 
exposed  animals.  For  each  of  the  behavioral  tests  (MWM,  EPM, 
and  OFT),  behavioral  outcomes  were  regressed  against  FA  val¬ 
ues  on  a  voxel-by- voxel  basis  to  identify  voxels  with  a  significant 
E-statistic  (p  <  0.05). 

A  complementary  region  of  interest  (ROI)  analysis  was  per¬ 
formed  to  assess  distributed  FA  abnormalities.  The  mean  FA  map 
from  all  animals  (n  =  60)  was  segmented  into  gray  and  white  mat¬ 
ter  using  a  threshold  of  0.25.  The  ROIs  were  manually  edited  to 
include  only  the  cerebral  white  matter  tracts  and  cerebral  cortex 
by  excluding  the  cerebellum,  brainstem,  and  subcortical  structures, 
and  each  hemisphere  was  assessed  separately.  The  FA  map  from 
each  animal  was  converted  to  z-scores  (z  =  {x—  (x)/o)  on  a  voxel- 
by-voxel  basis  where  x  is  the  individual  FA  value,  p,  and  a  are  the 
mean  and  standard  deviation  FA  value  of  the  sham  group,  respec¬ 
tively.  To  ensure  robustness,  an  additional  six  blast  sham  animals 
were  included  ( n  =  14).  Moreover,  a  leave-one-out  approach  was 
used  for  the  sham  group  by  calculating  the  z-scores  for  each  of  the 
individual  brains  compared  to  the  mean  and  standard  deviation  of 
the  remaining  13.  The  percentage  of  voxels  exceeding  a  z-score  of 
±2  was  calculated  separately  for  the  WM  and  GM  ROIs  (31,  34). 
The  resulting  quantitative  metric,  percentage  of  abnormal  FA  vox¬ 
els,  was  compared  to  the  sham  group  for  each  of  the  blast-exposed 
groups  using  a  Student’s  f-test  and  a  Bonferroni  correction  for 
multiple  comparisons. 
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HISTOLOGY 

Fixed  brains  were  processed  for  routine  paraffin  embedding  and 
cut  at  a  thickness  of  4  pm.  Selected  sections  were  stained  in  0.01% 
toluidine  blue  O  for  5  min  and  mounted.  Alternatively,  sections 
for  immunofluorescent  staining  were  placed  in  10  mM  sodium 
citrate  buffer  in  PBS,  microwaved  for  5  min,  and  allowed  to  cool 
to  room  temperature.  Sections  were  incubated  at  4°C  overnight 
with  primary  antibodies  for  glial  fibrillary  acidic  protein  (GFAP; 
1:500;  EMD  Millipore,  Billerica,  MA,  USA)  or  cleaved  caspase- 
3  (1:500,  EMD  Millipore,  Billerica,  MA,  USA).  Incubation  with 
the  secondary  antibody  (Alexa  Fluor;  1:500,  Molecular  Probes) 
was  performed  for  30  min  at  room  temperature.  Sections  were 
mounted  and  imaged  with  a  Nikon  microscope. 

RESULTS 

BLAST  SHOCKWAVE  AND  TBI  MODEL  CHARACTERIZATION 

Pressure  time  histories  depict  the  standard  experimental  approx¬ 
imation  of  the  Friedlander  waveform  of  the  primary  shockwave 
overpressure  (Figure  IB).  Peak  reflected  overpressures  measured 
at  the  rat  head  had  a  high  accuracy  and  reproducibility,  with  mea¬ 
sured  values  of  451  ±25  (n  =  21)  and  112±  14  (n  =  22)  kPa  for 
the  450  and  100 kPa  conditions,  respectively  (Figure  1C).  Posi¬ 
tive  durations  were  0.34  ±0.15  ms  for  the  high-blast  exposures 
and  0.46  ±0.21  ms  for  the  low-blast  exposures.  Pressures  mea¬ 
sured  near  the  rat  thorax  were  considerably  lower  at  63  ±  26  and 
19  ±  5  kPa,  respectively.  To  ensure  the  presence  of  the  rat  did 
not  alter  blast  pressure  measurements,  separate  blasts  were  per¬ 
formed  with  sensors  placed  in  the  location  of  the  rat  head  but 
without  a  rat  in  position.  Measurements  demonstrated  similar 
reflected  pressures  of  450  ±  23  kPa  (n  =  3)  and  incident  pressures 
of  1 10  ±  5.6  kPa  (n  =  6).  The  experimental  apparatus  also  limited 
head  motion.  Using  high-speed  video  analysis,  a  maximum  lateral 
displacement  of  6.98  mm  was  measured  in  the  high-blast  condi¬ 
tion,  with  an  angular  acceleration  of  <150krad/s2  (Figure  2B). 
Minimal  head  motion  was  noted  for  the  low-blast  condition  with 
an  acceleration  of  <10krad/s2.  Therefore,  any  behavioral  deficits 
or  structural  brain  damage  resulting  from  these  exposures  was 
most  likely  the  result  of  shockwave  passage  through  the  brain 
tissues  and  its  subsequent  effects  rather  than  head  rotational 
acceleration. 

MEMORY  DEFICITS  IN  SHOCKWAVE  EXPOSURE  GROUPS 

On  the  MWM  Visuo-Spatial  Learning  Paradigm,  all  experi¬ 
mental  groups  demonstrated  the  ability  to  learn  (Figure  3A), 
since  latency  to  find  a  hidden  platform  decreased  over  time 
(p  <  0.001).  On  the  first  day  of  testing,  there  were  no  sig¬ 
nificant  differences  in  latency  between  the  groups  (p  =  0.18). 
On  the  second  day,  a  main  effect  of  blast  severity  approached, 
but  did  not  meet  the  threshold  for  significance  (p  =  0.06).  On 
the  third  day,  the  main  effect  of  blast  severity  was  signifi¬ 
cant  (p  =  0.034).  Both  the  low-blast  (p  =  0.015)  and  high-blast 
(p  =  0.016)  groups  had  significantly  greater  latencies  to  find  the 
platform  than  the  sham  group  (Figure  3B),  but  the  two  groups 
were  not  significantly  different  than  one  another  ( p  =  0.97).  Time 
after  blast  had  no  effect  on  latency  between  the  4  and  30  DPI 
groups  {p  =  0.86),  indicating  persistent  cognitive  deficits  following 
shockwave  exposure. 


Session  1  Session  2  Session  3 


□  Blast  Sham  □  Low-blast,  30  DPI  ■  High-blast,  30  DPI 

□  Low-blast,  4  DPI  □  High-blast,  4  DPI 
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FIGURE  3  |  Morris  water  maze.  Memory  and  learning  deficits  were 
evident  in  the  blast-exposed  animals  compared  to  sham  animals  (A).  Swim 
path  traces  from  the  final  trial  (B),  with  each  animal  shown  in  a  different 
color,  reveal  the  difficulties  in  finding  the  hidden  platform  (gray)  in  the 
blast-exposed  animals  compared  to  the  sham  animals.  The  entry  point  is 
indicated  with  an  arrow.  Error  bars  indicate  SEM.  *p  <  0.05. 


The  number  of  unsuccessful  trials  across  all  18  sessions  of  the 
MWM  also  indicated  impaired  memory  and  cognition  in  animals 
exposed  to  blast  (Figure  3A).  There  was  no  significant  interac¬ 
tion  between  blast  severity  and  time  {p  =  0.73),  but  a  significant 
main  effect  of  blast  severity  was  evident  {p  <  0.019).  The  high- 
blast  group  had  significantly  more  unsuccessful  trials  compared 
to  the  sham  blast  group  (p  =  0.006),  whereas  the  low-blast  group 
was  not  significantly  different  from  animals  exposed  to  either  high- 
blast  (p  =  0.075)  or  sham  (p  =  0.15)  blast.  The  main  effect  of  time 
after  injury  for  the  number  of  unsuccessful  trials  was  not  signifi¬ 
cant  (p  =  0.18).  Time  spent  in  the  target  quadrant  had  a  significant 
interaction  between  blast  severity  and  time  after  injury  (p  =  0.04). 
In  the  post  hoc  analysis,  the  low-blast  group  at  30  DPI  spent  sig¬ 
nificantly  more  time  in  the  target  quadrant  than  the  high-blast 
groups  at  either  4  (p  =  0.029)  or  30  (p  =  0.01)  DPI.  No  other 
individual  groups  were  significantly  different  from  one  another.  A 
main  effect  of  blast  severity  was  also  significant  (p  =  0.0 1 ) ,  with  the 
low-blast  group  spending  more  time  in  the  target  quadrant  than 
the  high-blast  (p  =  0.008)  or  sham  (p  =  0.031)  groups.  All  other 
comparisons  for  time  in  the  target  quadrant  were  not  significant. 

OPEN  FIELD  ASSESSMENT 

The  OFT  was  used  to  assess  locomotor  activity  and  anxiety  behav¬ 
ior  ( Figure  4A ) .  Activity,  as  measured  by  the  total  distance  traveled 
during  the  test,  did  not  have  a  significant  interaction  between 
blast  severity  and  time  after  injury  (p  =  0.69).  A  significant  main 
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FIGURE  4  |  Open  field  test  and  elevated  plus  maze  Emotionality  and 
anxiety  was  examined  with  the  Open  Field  Test  (A)  and  the  Elevated  Plus 
Maze  (B)  for  each  of  the  four  blast  conditions  and  sham  group.  Error  bars 
indicate  SEM.  *p  <  0.05. 


effect  of  blast  severity  was  observed  ( p  =  0.029).  Compared  to  the 
sham  blast  group,  animals  exposed  to  low  (p  =  0.018)  or  high 
(p  =  0.0 1 1 )  blast  conditions  had  significantly  less  activity,  but  they 
were  not  significantly  different  from  one  another  {p  =  0.76).  The 
effect  of  time  after  injury  was  not  significant  (p  =  0.84).  Likewise, 
the  amount  of  time  spent  in  the  center  of  the  OFT,  a  measure  of 
anxiety,  did  not  have  a  significant  interaction  between  blast  sever¬ 
ity  and  time  after  injury  {p  =  0.20).  Time  in  the  center  was  also  not 
significant  for  either  the  main  effect  of  blast  severity  (p  =  0.51)  or 
time  after  injury  {p  =  0.63). 

EMOTIONALITY  CHANGES  IN  SHOCKWAVE  EXPOSURE  GROUPS 

Exploratory  behavior  on  the  EPM  served  as  an  indicator  of  emo¬ 
tionality.  The  amount  of  time  spent  in  the  open  arms  of  the  EPM 
(Figure  4B),  a  metric  associated  with  anxiety,  did  not  have  a  sig¬ 
nificant  interaction  with  blast  severity  and  timing  after  injury 
{p  =  0.39).  Elowever,  the  main  effect  of  blast  severity  was  signifi¬ 
cant  {p  =  0.037).  The  high-blast  group  spent  significantly  less  time 
in  the  open  arms  compared  to  the  low-blast  ( p  =  0.013)  group.  The 
high-blast  (p  =  0. 1 1 )  and  low-blast  (p  =  0.71 )  groups  were  not  sig¬ 
nificantly  different  from  the  sham  group.  The  main  effect  of  time 
after  injury  was  not  significant  (p  =  0. 19). 

The  number  of  arm  changes  in  the  EPM  did  not  have  a  sig¬ 
nificant  interaction  between  blast  severity  and  time  after  injury 


(p  =  0.63).  The  main  effect  of  blast  severity  was  not  significant 
(p  =  0.35),  although  the  number  of  arm  changes  was  decreased 
by  26  and  29%  in  the  low  and  high-blast  groups,  respectively 
(Figure  4B).  The  main  effect  of  timing  after  injury  was  not 
significant  (p  =  0.98). 

SHOCKWAVE  EXPOSURE  CAUSED  SPATIALLY  DEPENDENT  AND 
PROGRESSIVE  DTI  ABNORMALITIES 

Shockwave  exposure  caused  significant  changes  in  FA.  A  distinct 
pattern  of  abnormalities  emerged  related  to  blast  severity  and  time 
after  injury,  with  significant  differences  between  each  of  the  blast- 
exposed  groups  compared  to  the  sham  group  (Figure  5A).  At  4 
DPI,  significant  decreases  in  FA  were  evident  in  the  ipsilateral  cor¬ 
tex,  medial  prefrontal  cortex  (mPFC),  hippocampus,  and  thalamus 
and  were  more  pronounced  in  the  animals  exposed  to  high-blast 
than  those  exposed  to  low-blast  shockwaves.  At  30  DPI,  significant 
decreases  in  FA  were  evident  in  both  blast  exposures.  Brain  regions 
affected  at  30  DPI  included  those  affected  at  4  DPI  along  with 
additional  injury  evident  in  the  contralateral  cortex  and  brain¬ 
stem.  The  voxelwise  statistical  tests  were  uncorrected  for  multiple 
comparisons  to  aid  in  visualizing  the  spatial  pattern  of  FA  changes 
across  all  of  the  blast-exposed  groups  with  the  same  sham  group 
used  for  all  comparisons.  The  decreases  in  FA  were  greater  and 
more  widespread  in  the  group  exposed  to  high-blast  compared 
to  those  exposed  to  low-blast.  The  three-dimensional  visualiza¬ 
tion  of  significant  differences  are  shown  in  Figure  5B.  The  effect 
sizes  (Cohen’s  D)  for  the  significant  clusters  in  each  of  the  groups 
compared  to  the  sham  group  were  0.417  and  0.407  for  the  low- 
blast  and  high-blast  at  4  DPI,  respectively,  and  0.493  and  0.495 
for  the  low-blast  and  high-blast  groups,  respectively,  at  30  DPI. 
When  corrected  for  multiple  comparisons  at  p  <  0.05  by  control¬ 
ling  the  familywise  error,  the  voxelwise  differences  did  not  survive 
and  were  therefore  followed  up  by  an  unbiased  ROI  analysis. 

A  ROI  analysis  was  performed  to  complement  the  voxelwise 
findings  and  quantitatively  examine  the  laterality  of  the  DTI 
changes  (Figure  6).  ROIs  encompassing  the  left  (ipsilateral)  or 
right  (contralateral)  cerebral  white  matter  (Figure  6A)  or  cortical 
gray  matter  (Figure  6B)  revealed  that  mean  FA  was  not  signifi¬ 
cantly  different  between  any  of  the  groups.  However,  an  alternative 
method  quantifying  the  number  of  abnormal  FA  voxels  for  each 
animal  relative  to  a  spatially  aligned  sham  group  (z-score)  revealed 
significant  differences  between  groups.  In  the  white  matter,  the 
percentage  of  abnormal  (Izl  >  2)  voxels  was  significantly  greater  in 
all  of  the  blast-exposed  groups  compared  to  shams  in  the  ipsilateral 
white  matter,  with  greater  changes  observed  at  30  DPI  (Figure  6C). 
In  the  contralateral  white  matter,  no  significant  differences  were 
observed.  However,  in  the  cerebral  cortex  (Figure  6D),  the  number 
of  abnormal  FA  voxels  was  significantly  greater  in  the  ipsilateral 
cortex  for  all  blast-exposed  groups  compared  to  the  sham  group 
(all  p  <  0.05,  corrected  for  multiple  comparisons).  In  the  con¬ 
tralateral  hemisphere,  the  low-blast  (p  =  0.0009)  and  high-blast 
(p  =  0.0005)  exposure  groups  at  the  chronic  time  points  were  sig¬ 
nificantly  different  than  the  sham  group,  whereas  no  significant 
different  differences  were  noted  at  the  acute  time  point.  In  general, 
the  ROI  results  confirmed  the  previous  voxelwise  observations.  To 
identify  whether  the  FA  changes  were  predominantly  associated 
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FIGURE  5  |  Microstructural  abnormalities  following  primary  blast  The 

brains  of  animals  exposed  to  blast  shockwave  underwent  DTI  and  compared 
to  animals  exposed  to  sham  blast  (A).  Significant  FA  decreases  were  evident 
at  4  (left)  and  30 days  post-blast  (right).  At  both  timepoints,  greater  numbers 
of  significant  voxels  were  evident  in  the  animals  exposed  to  high-blast 
compared  to  those  at  exposed  to  low-blast.  Compared  to  4  days  post-blast, 


there  were  more  significant  voxels  at  30  days  post-blast  in  both  blast 
exposures.  Decreased  FA  was  most  prominent  in  the  cortex,  but  also 
included  the  thalamus,  and  ipsilateral  ventral  hippocampus.  The 
three-dimensional  distribution  of  significant  changes  compared  to  sham  are 
shown  (B)  as  viewed  from  the  top  of  the  brain.  The  direction  of  blast  is 
indicated  by  the  arrow  and  is  identical  for  both  panels. 


with  changes  in  axial  diffusivity  or  radial  diffusivity,  these  values 
were  obtained  from  the  abnormal  FA  voxels  from  each  animals 
and  averaged  across  groups.  There  were  no  significant  differences 
in  AD  or  RD  between  any  of  the  groups  compared  to  the  sham 
group. 

RELATIONSHIP  BETWEEN  FRACTIONAL  ANISOTROPY  AND  BEHAVIOR 

Linear  regression  analysis  was  performed  to  identify  relationships 
between  behavioral  outcomes  and  microstructure  assessed  with 
DTI  (Figure  7).  In  the  MWM,  clusters  of  voxels  with  significant 
linear  relationships  between  FA  and  behavioral  outcomes  were 
evident  bilaterally  in  the  hippocampus,  predominantly  in  the  dor¬ 
sal  region.  A  significant  region  was  also  evident  in  the  motor 
cortex,  predominantly  in  the  contralateral  cortex.  Across  all  sig¬ 
nificant  voxels,  target  quadrant  percentage  was  a  better  predictor 
of  FA  (p  <  0.001)  than  number  of  unsuccessful  trials  ( p  =  0.045) 
or  latency  (p  =  0.088).  In  the  OFT,  outcomes  were  significantly 
related  to  FA  of  the  contralateral,  ventral  hippocampus,  and  a  small 
portion  of  the  ipsilateral  external  capsule.  FA  in  these  regions  was 
significantly  associated  with  both  outcomes  of  the  OFT,  time  in  the 
center  ( p  <  0.001),  and  distance  traveled  (p  <  0.001).  In  the  EPM, 
only  a  small  region  of  significance  was  located  in  the  contralateral 
cortex.  Neither  the  number  of  arm  changes  (p  =  0.082)  or  time 
spent  in  the  open  arms  ( p  =  0.30)  were  independent  predictors  of 
FA  in  this  region. 

HISTOLOGICAL  EVIDENCE  OF  INJURY 

No  hemorrhage  was  evident  in  any  of  the  animals  exposed  to 
shockwaves  based  on  visual  inspection  of  the  brain  surface.  A 


qualitative  analysis  revealed  greater  astrocyte  hypertrophy  in  the 
mPFC  (Figure  8),  a  region  exhibiting  significant  DTI  abnor¬ 
malities  (Figure  4A),  following  high-blast  shockwave  exposure 
compared  to  low-blast  shockwave  exposure.  A  similarly  greater 
hypertrophy  was  evident  at  chronic  compared  to  acute  time 
points.  Correspondingly,  sections  stained  with  Toluidine  blue 
revealed  chromalytic  neurons  in  the  cortex  of  brains  exposed  to 
high-blast  shockwaves  at  both  acute  and  chronic  time  points, 
whereas  the  brain  exposed  to  low-blast  shockwaves  contained 
only  sparse  darkly  stained  neurons.  Staining  for  cleaved  caspase- 
3,  a  nuclear  marker  of  apoptosis,  revealed  numerous  apop- 
totic  cells  within  the  ipsilateral  cortex,  whereas  no  apoptotic 
nuclei  were  evident  in  the  contralateral  cortex  or  at  30  days 
post-blast. 

DISCUSSION 

Animal  brain  injury  models  have  proven  essential  to  understand 
the  consequences  of  blast  injury.  Although  the  precise  mechanisms 
of  blast  neurotrauma  are  the  subject  of  considerable  investigation 
and  debate,  several  important  experimental  factors  may  have  a 
considerable  influence  on  the  findings.  The  current  investigation 
employed  an  experimental  rodent  model  designed  to  be  a  trans¬ 
lational  paradigm  representative  of  human  blast  injury  dynamics. 
In  our  animal  model,  the  body  of  each  rat  was  protected,  so  only 
the  head  was  exposed  to  shockwaves.  This  protective  approach  is 
consistent  with  the  setting  of  ballistic  personal  protective  body 
armor  worn  by  military  personnel  being  effective  at  attenuat¬ 
ing  blast  overpressure  exposure  to  the  thorax  (40).  Field  injuries 
are  more  likely  attributable  to  shockwave  exposure  of  the  head, 
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FIGURE  6  |  Region  of  interest  DTI  analysis  The  mean  FA  of  the  white 
matter  tracts  (A)  or  cerebral  cortex  (B),  assessed  separately  for  the 
ipsilateral  and  contralateral  hemispheres,  were  not  significantly 
different  between  any  of  the  groups.  In  contrast,  the  percentage  of 
voxels  with  FA  values  >2  standard  deviations  from  the  mean  (z-scores) 
compared  to  the  sham  group  revealed  significant  group  differences.  In 
the  white  matter  (C)  no  significant  differences  were  evident  in  either 


the  ipsilateral  or  contralateral  hemispheres.  In  the  cerebral  cortex 
(D),  the  percentage  of  abnormal  voxels  ipsilateral  to  the  blast  exposure 
was  significantly  greater  in  all  blast-exposed  groups  compared  to  the 
sham  group.  In  the  contralateral  hemisphere,  only  groups  30 days 
post-blast  had  a  significantly  greater  percentage  of  abnormal  FA  voxel 
compared  to  the  sham  group.  Error  bars  indicate  standard 
deviation. *p  <  0.05  compared  to  sham. 


and  injury  tolerances  are  different  between  shockwave  exposures 
to  the  head  and  brain  tissue  compared  to  those  involving  pul¬ 
monary  injury.  The  scaling  of  shockwave  characteristics  from  the 
rat  to  the  human  is  another  important  experimental  concern. 
Construction  of  explosive  devices  using  artillery  shells  (41)  has 
demonstrated  overpressure  durations  typically  <10  ms.  Scaling 
ratios  presented  by  Bass  et  al.  (42),  albeit  for  pulmonary  expo¬ 
sures,  can  be  used  to  relate  overpressure  durations  from  rodent 
models  to  human-equivalent  values.  In  the  current  study,  scaled 
overpressure  durations  were  consistent  with  human  exposures  of 
2-3  ms.  Isolating  the  effects  of  blast  shockwave  from  head  accel¬ 
eration  is  an  important  concept  to  understand  the  mechanisms 
of  blast  injury.  In  the  current  study,  the  head  was  laterally  sup¬ 
ported  during  shockwave  exposure,  which  prevented  appreciable 
head  excursion  during  the  low-blast  exposures  and  limited  excur¬ 
sions  during  high-blast  exposures.  Angular  accelerations  were  well 
below  the  threshold  for  mild  concussion  in  the  rat  (43).  Similarly, 
positioning  rats  outside  and  off-axis  from  the  shock  tube  limited 
exposure  to  venting  gas  (i.e.,  blast  wind).  While  exposure  outside 
of  the  tube  introduces  shockwave  dynamics  such  as  diffraction 
and  weakening,  all  rats  were  exposed  to  a  repeatable  shockwave 
profile  (Figure  1 )  without  blast  wind  confirmed  by  measured  pres¬ 
sure  profiles.  Collectively,  the  resulting  behavioral,  imaging,  and 
histological  consequences  of  blast  injury  were  attributable  to  the 
primary  shockwave  exposure  and  not  pulmonary  or  rotational 
acceleration  mechanisms. 


Primary  blast  injury  caused  significant  memory  deficits  in 
shockwave-exposed  rats  compared  to  sham-exposed  animals,  with 
more  pronounced  deficits  in  the  higher  shockwave  exposure 
group.  These  data  provide  preliminary  support  for  the  theoret¬ 
ical  assumption  of  a  “dose  effect”  associated  with  more  severe 
or  repetitive  blast  TBI  exposure.  The  observed  memory  and 
cognitive  deficits  were  generally  consistent  with  findings  from 
other  experimental  investigations  of  shockwave-induced  TBI  in 
rodents  (21,  44),  although  the  magnitude  of  blast  exposure  dif¬ 
ferences  between  the  studies  may  have  implications  regarding 
injury  tolerance  to  blast.  Goldstein  and  colleagues  (14)  demon¬ 
strated  significant  cognitive  deficits  after  77  kPa  exposures,  but  the 
high  magnitude  head  rotational  accelerations  potentially  caused 
injury  through  various  mechanisms.  The  cognitive  deficits  in 
the  present  protocol,  which  incorporated  head  restraint,  were 
directly  attributable  to  the  interaction  of  the  shockwave  and 
the  head.  A  notable  difference  between  the  present  study  and 
previous  investigations  is  the  presence  of  significant  cognitive 
deficits  acutely  after  injury.  Increased  latency  to  find  the  plat¬ 
form  and  greater  numbers  of  unsuccessful  trials  were  evident  for 
shockwave-exposed  groups  at  <4  days  post  injury.  Other  studies 
have  reported  no  differences  between  controls  and  shockwave- 
exposed  rats  at  3-8  days  post  injury  (21,  44,  45).  The  disparity 
could  be  related  to  the  use  of  the  Barnes  Maze  assessment  in 
other  studies,  compared  to  the  use  of  the  Morris  Water  Maze  in 
the  current  study,  since  the  two  assessments  may  have  different 
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FIGURE  7  |  Diffusion  tensor  imaging  correlations  with  behavior.  Across 
all  animals  exposed  to  blast,  multiple  linear  regression  was  used  to 
determine  whether  behavioral  outcomes  were  associated  with  brain 
microstructure.  For  each  behavioral  test,  the  outcome  measures  obtained 
from  it  were  regressed  against  the  FA  value  from  all  animals  using  a  F-test, 
controlling  for  blast  severity  and  DPI.  In  the  MWM,  FA  was  significantly 
related  to  behavioral  performance  in  the  hippocampus  and  portions  of  the 
motor  cortices.  In  the  OFT,  a  region  of  the  ventral  hippocampus  and  region 
of  the  external  capsule  were  significantly  related  to  the  behavioral 
performance.  The  EPM  was  related  to  FA  in  only  a  small  region  of  the 
contralateral  cortex. 


sensitivities.  These  differences  may  also  be  attributable  to  differ¬ 
ences  in  the  injury  protocol  such  as  the  used  of  head  restraints  and 
torso  protection. 

Neurobehavioral  and  emotional  changes  assessed  using  the 
EPM  and  OFT  were  also  generally  consistent  with  previous  reports 
(44).  Elevated  levels  of  anxiety  were  evident  for  shockwave- 
exposed  rats,  with  some  indication  of  a  dose  response  and  transient 
symptoms.  Sham  rats  demonstrated  the  greatest  magnitude  of 
open  arm  time,  with  the  groups  exposed  to  blast  spending  progres¬ 
sively  less  time  in  the  open  arms  at  the  acute  time  point.  Likewise, 
open  field  center  time  decreased  in  the  sham  group  compared  to 
the  groups  exposed  to  blast.  These  findings  indicate  higher  levels 
of  anxiety  in  more  severely  injured  rats,  again  supporting  a  dose 
effect  of  blast  TBI.  This  type  of  dose  dependency  for  emotional 
changes  following  shockwave  injury  may  have  clinical  implica¬ 
tions  for  the  understanding  of  patient  condition  and  outcomes. 
Whether  emotional  changes  are  transient  or  persistent  follow¬ 
ing  shockwave  exposure  has  been  mixed.  In  the  current  study, 
the  time  in  the  center  of  the  open  field  was  decreased  in  the 
high-blast  group  acutely,  but  recovered  by  1  month  post  injury. 
However,  open  arm  time  in  EPM  decreased  acutely  in  the  high- 
blast  group  and  did  not  resolve  by  1  month  post  injury.  Similar 
behavioral  paradigms  at  other  laboratories  have  observed  pro¬ 
gressive  emotional  changes  at  44  DPI  (20),  but  no  deficits  at 
66  DPI.  Conversely,  persistent  decreases  in  open  arm  time  from 


48  h  to  1  month  post  injury  have  also  been  observed  (44).  Emo¬ 
tional  deficits  have  also  been  shown  to  persist  up  to  6  months 
post- injury  (46).  The  disparities  maybe  related  to  the  interaction 
between  stress  and  emotional  dysfunction  (44,  47)  and  require 
further  investigation,  particularly  in  relation  to  the  development 
of  PTSD  in  military  personnel  and  returning  veterans  exposed  to 
blast  during  combat  (2,  4,  7). 

The  finding  of  DTI  abnormalities  and  their  spatial  distribu¬ 
tion  provides  some  insight  into  the  mechanisms  of  injury  fol¬ 
lowing  shockwave  exposure.  At  acute  time  points,  decreased  FA 
in  the  superficial  cortex  ipsilateral  to  the  shockwave  suggests  a 
local  effect  of  the  shockwave  on  the  brain  tissue,  and  is  consis¬ 
tent  with  previous  histological  studies  demonstrating  gliosis  in 
the  lateral  (14)  and  medial  prefrontal  cortices  (44).  The  effects 
may  be  related  to  high  shear  strain  or  microcavitation  at  the 
interface  between  the  brain  and  cerebrospinal  fluid  shown  in 
simulations  of  blasts  in  the  human  (48)  and  rat  brain  (27).  In 
the  current  study,  injuries  progressed  from  acute  to  chronic  time 
points  and  involved  greater  regions  of  the  cortex  as  well  as  por¬ 
tions  of  the  brainstem,  suggesting  an  ongoing  response  similar 
to  that  observed  elsewhere  (49).  The  injury  in  the  superficial 
cortex  could  potentially  be  related  to  the  acute  disruption  of 
the  blood-brain  barrier  (17)  and  it  sequelae.  Indeed,  gliosis  is  a 
prominent  pathology  after  blast  (15)  and  hypertrophic  astrocytes 
were  evident  in  the  mPFC  in  the  current  study  (Figure  7).  White 
matter  injury  was  also  evident  in  the  current  study  (Figure  4), 
consistent  with  previous  histological  reports  of  axonal  injury  in 
the  brainstem  and  cerebellum  (50).  Furthermore,  the  regions  of 
injury  could  be  associated  with  auditory  regions  in  the  brain 
(51).  Further  studies  will  be  necessary  to  understand  the  com¬ 
plex  relationship  between  biomechanical  and  other  consequences 
of  blast  injury. 

The  current  study  employed  voxelwise  analysis  of  FA  images 
since  the  orientation  of  the  blast  was  consistent  across  all  ani¬ 
mals.  However,  since  blast  neurotrauma  in  humans  is  highly 
variable  with  respect  to  intensity,  direction,  and  other  factors, 
alternative  approaches  to  capture  the  heterogeneity  have  been 
proposed.  Comparing  an  individual  to  a  control  group  to  derive 
group-averaged  measures  (31,  34)  or  single  subject  maps  (52)  of 
abnormal  FA  was  shown  to  have  greater  sensitivity  to  TBI  than 
measures  of  mean  FA.  In  the  current  study,  this  approach  was  used 
to  quantify  the  voxelwise  changes.  As  expected,  there  were  no  dif¬ 
ferences  in  whole  brain  mean  FA  values  between  any  of  the  groups. 
However,  quantifying  the  number  of  abnormal  FA  voxels  in  each 
subject  relative  to  the  control  group  and  group-averaging  this  met¬ 
ric  was  sensitive  to  group  differences  (Figure  6).  The  results  were 
largely  consistent  with  the  voxelwise  findings,  but  also  provide 
greater  insight  into  the  directionality  of  the  FA  changes,  since  this 
is  conflicting  in  human  DTI  studies  following  TBI.  In  the  cortex, 
the  FA  decrease  was  greater  with  higher  blast  intensity  and  with 
longer  durations  from  the  injury.  In  contrast,  FA  increased  in  the 
white  matter  acutely  and  appeared  to  normalize  at  the  30  day  time- 
point.  This  is  the  opposite  that  of  the  white  matter  FA  decreases 
observed  in  severe  TBI  (53).  In  human  mTBI,  however,  increased 
white  matter  FA  was  often  observed  in  the  acute  (54)  and  semi¬ 
acute  (34)  periods  whereas  decreased  FA  was  consistently  observed 
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FIGURE  8  |  Histological  evidence  of  blast  effects  Brain  sections  from  the 
regions  exhibiting  DTI  changes,  approximately  +2.0  mm  Bregma,  were 
stained  with  toluidine  blue  (A)  Sections  revealed  numerous  "dark"  neurons 
in  the  brains  exposed  to  high-blast  (C)  at  both  4  and  30  days  post-blast, 
particularly  in  the  ipsilateral  cortex,  but  fewer  in  the  contralateral  cortex.  Only 
sporadic  abnormal  neurons  were  observed  in  the  low-blast  and  sham  brains. 
In  a  brain  exposed  to  high-blast,  staining  for  cleaved  caspase-3  (B)  was 


observed  in  the  ipsilateral  cortex  at  4  days  post-blast,  but  not  in  the 
contralateral  cortex,  and  the  stain  was  no  longer  present  at  30  days  post-blast. 
No  staining  was  present  in  the  sham  brain.  In  the  medial  prefrontal  cortex,  a 
region  with  significant  FA  changes,  astrocytes  stained  with  GFAP  exhibited 
extensive  hypertrophy  in  the  high-blast  brain  at  30  days  compared  to  the 
sham  brain,  but  astrogliosis  was  noticeably  less  at  4 days  post-blast  or  in  the 
low-blast  condition. 


in  the  more  chronic  phases  (30-32).  Thus,  although  FA  is  a  sen¬ 
sitive  biomarker  of  injury  following  TBI,  interpreting  the  changes 
in  FA  are  complex  with  respect  to  the  underlying  pathology  (55). 

The  relationship  between  DTI,  cognitive,  and  neurobehavioral 
abnormalities  identifies  a  putative  link  between  regional  brain 
abnormalities  and  the  resulting  clinical  manifestations  of  blast 
injury.  While  the  current  findings  were  observed  in  an  experi¬ 
mental  animal  model,  they  likely  have  significant  translational 
relevance  to  our  understanding  of  the  pathophysiology  and  clinical 
effects  of  blast  injury  in  human.  Impairment  in  memory  function 
was  related  to  abnormalities  in  the  hippocampus  (Figure  7),  as 
expected  although  the  precise  mechanism  for  the  changes  is  not 
clear.  A  quantitative  assessment  of  the  biological  basis  of  the  DTI 
findings  would  provide  greater  insight  into  the  injuries  evident  in 
humans  after  blast.  One  limitation  of  the  current  study  is  the  use 
of  fixed  tissues.  Fixation  has  been  shown  to  maintain  anisotropy 
(56),  but  in  vivo  measurements  would  enable  a  better  appreci¬ 
ation  of  longitudinal  changes  following  TBI.  On  the  other  hand, 
imaging  of  fixed  tissues  affords  substantial  improvements  in  image 
quality  and  resolution  compared  to  in  vivo  and  thereby  better  sen¬ 
sitivity  to  subtle  changes.  Another  limitation  is  that  blast  injuries 
were  performed  at  a  single  orientation  relative  to  the  shock  tube 
(side-on).  The  importance  of  orientation  has  been  demonstrated 
using  simulations  (27,48)  and  experimental  studies  (19,57),buta 
whole  brain  assessment  such  as  that  provided  by  DTI  has  yet  to  be 


performed.  Utilization  of  serial,  in  vivo  MRI  will  be  important  to 
fully  characterize  the  effects  of  blast  TBI  and  the  evolution  from 
acute  to  chronic  injury. 

CONCLUSION 

In  a  rat  model  of  blast  TBI,  primary  blast  exposure  was  sufficient  to 
induce  cognitive  and  neurobehavioral  deficits  and  microstructural 
injury  detected  with  DTI.  Importantly,  both  outcomes  scaled  with 
the  magnitude  of  the  blast  shockwave,  reflecting  a  dose  effect  that 
has  implications  to  the  fundamental  pathophysiology  and  neuro¬ 
logic  sequelae  of  blast  TBI  in  humans.  Moreover,  the  correlation 
between  neurological  deficits  and  DTI  findings  underscores  the 
connection  between  brain  injury  and  behavior  that  has  particular 
relevance  to  advancing  the  science  of  therapeutic  interventions  and 
prevention  of  long-term  effects  of  blast  injury  in  military  service 
members  and  veterans. 
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